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ABSTRACT 



Aims. We develop a novel way of finding total mass density profiles in Sersic ellipticals, to about 3 times 
the major axis effective radius, using no other information other than what is typically available for distant 
galaxies, namely the observed surface brightness distribution and the central velocity dispersion ctq. 
Methods. The luminosity density profile of the observed galaxy is extracted by deprojecting the measured 
brightness distribution and scaling it by a fiduciary, step-function shaped, raw mass-to-light ratio profile 
(MIL). The resulting raw, discontinuous, total, 3-D mass density profile is then smoothed according to a 
proposed smoothing prescription. The parameters of this raw M/L are characterised by implementing the 
observables in a model-based study. 

Results. The complete characterisation of the formalism is provided as a function of the measurements of 
the brightness distribution and ctq. The formalism, thus specified, is demonstrated to yield the mass density 
profiles of a suite of test galaxies and is successfully applied to extract the gravitational mass distribution in 
NGC 3379 and NGC 4499, out to about 3 efi'ective radii. 

Key words. Methods: analytical - Galaxies: fundamental parameters (masses) 



1. Introduction 

Any evaluation of the total mass in distant galaxies 
is a struggle against the paucity of available obser- 
vational evidence. Photometry is hardly enough to 
indicate the content of both the luminous as well 
as the dark matter, unless the functional dependence 
between luminosity content and dark mass is ac- 
cessible. This is of course not the case; the exis- 
tence of such a relation is itself uncertain, espe- 
cially in early type galaxies. While kinematic in- 
formation of tracers has often been advanced as in- 
dicators of mass distributions in galaxies, the im- 
plementation of such information is tricky, pri- 
marily because of the mass-anisotropy degeneracy. 
Thus, one often resorts to cleverly designed ob- 
servational techniques and/or algorithms and for- 
malisms which thrive even in light of the lim- 
ited measurements. Examples o f these include the 
Plane tary Nebula Spectrograph ("Romanowskv et qZj 
l2004l) . the NM AGIC code (de Lorenzi ef a/. ^OOTb 
and CHASSIS dChakrabartv & Sahall2001l) . 



The X-ray emission from X-ray active systems 
can be analysed to offer insight into the gravitational 
mass distribution in the galax y, under t he assumption 
of hydrostatic eqttilibrium ({ Humphrey et al. 2008; 



Lemze et al.1l2008HMahdavi et al. 2008; Zhang etal 
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2007| iFukazawa et aTl 120061; oiSullivan & Ponma^ 
2004, to cite some recent work). However, what 
makes this method potentially unreliable is the lack 
of information about the distribution of the frac- 
tion of hot gas that is in hydrostatic equ ilibrium 
(IChurazov et al.ll2008t loiehl & Stadeiil2007h . 

Comparatively, a more stable route to mass dis- 
tribution determination is via lensing measurements. 
However, the biggest shortcoming of mass determi- 
nation from lensing measurements alone is the un- 
availability of the full three-dimensional mass dis- 
tribution. To improve upon this, lensing data is of- 
ten supple mented by dynami cally obtained mass 
estimates dCzoske et all 120081; fe olton et al] l2008t 
iGavazzi et al.ll2007l; iKoopmans & Treuu200f ). 

However, there are questionable implementa- 
tional problems involved in (parametric) dynamical 
mass determination, the chief of which are typically 
the mass-anisotropy degeneracy, binning-triggered 
instability of scant velocity dispersion data, reliance 
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on the modelling of the stellar mass density and 
an even more fundamental worry caused by the as- 
sumption of one smooth parametric representation 
of the phase-space distribution function of the used 
tracer and by its relation to the phase-space density 
of the whole galaxy. This is of course in addition 
to the uncertainties in the mass distribution resulting 
from substituting the real geometry of the system by 
sphericity, as is typically done with all mass determi- 
nation procedures. Above al l, using tracer kinemat- 
ics for mass determ ination jde Lorenzi et al. 
iDouglas et. all 120071 etc.) is limited in applicability 
given the reliance on the size of the tracer kinematic 
data! Large data sets are of course hard to attain in 
systems that are not close by. Moreover, this method 
is unsuitable for fast evaluation of the mass distri- 
bution of individual galaxies that are members of a 
large sample, as for example, a galaxy obseved in a 
large survey. 

On the contrary, it would be highly beneficial 
to design a method that is comparatively less data- 
intensive in that it demands only what is easily avail- 
able from observations. We advance a methodology 
that provides total gravitational mass density distri- 
butions to about 3 e ffective radii, as compa red to only 
1 effective radius (|Cappellari et al.l l2006h . in a fast 
and easy-to-implement fashion. 

This advanced f ormahsm is inspired by a trick 
that was reported in IChakrabartvl ( |2007|) (hereafter. 
Paper I). This trick involves the exploitation of only 
photometry and the central velocity dispersion mea- 
sure {(To) in a galaxy, in order to generate the to- 
tal local mass-to-light ratio (M/L) profile to a dis- 
tance that is about thrice the semi-major axis effec- 
tive radius. This cutoff distance is described in details 
below. The prescription for constructing this profile 
was provided in Paper I, though only for a certain 
class of power-law galaxies. However, the exact na- 
ture of this prescription is very much a function of the 
photometric class that the galaxy belongs to. Thus, 
the formula reported in Paper I cannot be invoked to 
shed light on the mass distribution in galaxies that 
betray a different (and more ubiquitous) photometric 
class, eg. ellipticals, the surfac e brig h tnes s of which 
can be fit by a Sersic profile (ISersiall968l) . This is 
precisely what is reported in this paper. 

This paper is arranged as follows: the basic 
framework of the suggested formalism is discussed 
in Section 2, followed by a note on the models 
that we use. The method used to obtain the sought 
functional forms is briefly mentioned in Section 4. 
Results obtained from our work are subsequently 
discussed in Section 5. Tests of the method are de- 
scribed in Section 6 while Section 7 deals with appli- 
cations to real galaxies NGC 3379 and NGC 4494. 
The paper is rounded off with a section devoted to 
discussions of relevant points. 



2. Formalism 

The only Sersic model that was considered in Paper I 
did actually indicate that the mass estimation trick 
suggested for the power-law systems might be pos- 
sible for Sersic galaxies too. Following this lead, as 
in Paper I, we first invoke a raw two-stepped M/L 
profile of the Sersic galaxy at hand, where it is the 
distribution of M/L with the major axis coordinate x, 
that is relevant. This raw M/L distribution is subse- 
quently smoothed, (according to the smoothing pre- 
scription provided in Paper I and discussed below) to 
provide the real M/L distribution of the system, to a 
distance that is by definition, 3 times the major-axis 
effective radii for a Sersic galaxy with sersic index 
n=4 but is an approximation for the major-axis effec- 
tive radius for all other values of n. The formal def- 
inition of this distance is given in Equation [1] while 
the justification for our choice of this length scale 
is delineated in Section 18.41 Figure [T] represents a 
schematic diagram of this raw M/L profile against 

X. 

As can be appreciated from Figure [1] the unpro- 
cessed M/L profile is a two-stepped function that is 
characterised by three free parameters, including the 
position of the step, or the jump radius Xi„ and the 
M/L amplitudes inside and outside x,„. The object of 
the current exercise is to fully characterise this un- 
processed M/L profile, so that upon smoothing, the 
final, total, local M/L distribution over x is retrieved. 
This final form of the M/L distribution is signifi- 
cantly different from the discontinuous, 2-stepped, 
raw M/L that we initially choose to work with. In 
fact, the final form of the profile is smooth and the 
M/ L values vary a lot with x, sometimes abruptly, de- 
pending on the details of the model galaxy. Figure 4 
in Paper I represents a comparison of one such raw 
M/L profile and the final smoothed M/L distribution 
that is advanced as the representative M/L profile for 
the test system at hand. In fact, later in Figure|6l a raw 
M/ L (exemplified in Figure[T]) is compared to the true 
M/L distribution of the test galaxy under considera- 
tion - the difference between the raw and final forms 
of the M/L profiles is clear in that figure. 

As in Paper I, x,,, is set equal to 3Xc, where 



and m is the slope of the straight line that is fit to 
the plot of core-removed logio(I) against x^^'^ (as 
in Paper I), i.e. is the major axis equivalent of 
the effective radius, if the Sersic index is 4. For all 
other values of the Sersic index, is at most an ap- 
proximation to the major axis effective radius. The 
advantage of using X^ over the exact definition of 
the major-axis effective radius is described below in 
SectionlSH 

The smoothing prescription used here is the same 
as in Paper I - we smooth the raw total mass density 
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tral value of 1000 Lopc ^. Thus, we want to find the 
function f{n, X^, T,„), where 



Fig. 1. Schematic view of a raw M/L distribution 
along the major axis coordinate x, for a model galaxy. 
The dashed line marks the position of the jump ra- 
dius Xi„ which as defined in Equation [T] is 3 times 
the major axis eff'ective radius of the model galaxy, 
if the galaxy is described by a Sersic index of 4; oth- 
erwise, x,„ is an approximation for the major-axis ef- 
fective radius. The amplitudes inside and outside X,„ 
have been marked as T,„ and Tom- While the mea- 
sured central velocity dispersion is used to estimate 
the range of T,„ values allowed, To„, is constrained 
via the sought analytical relation that connects Tout 
to T,„. The point of this paper is to seek such rela- 
tions, in order to fully characterise the raw M/L pro- 
file; once the raw M/L distribution is known, it is 
smoothed to obtain the final, total M/L ratio profile 
of the galaxy, which is significantly different from 
this unprocessed M/L shown above. 



profile by two successive applications of a box filter 
of size corresponding to X,,. 

The amplitude of M/L(x) for x < Xj,, is referred 
to as Yin and that for x > Xi„ it is Tom- Now that Xin 
has been pinned down by construction, we hope to 
get a constraint on the available choices for T,„ from 
the central velocity dispersion (ctq) and also hope 
to identify a functional dependence of To„i on Y,„. 
The exact form of such a constraint or function is 
yet unknown but we begin by expecting these to be 
defined in terms of the photometric parameters that 
describe the surface brightness profile of the Sersic 
galaxy under consideration, namely, the Sersic index 
n and X,, that we have defined above (Equation [U. 
The central brightness is not a free parameter since 
we normalise all luminosity density profiles to a cen- 



"^out — f{n,Xg, Tin) 



(2) 



The inspiration for the hypothesis that / depends 
on the photometric parameters is discussed below. 

Given that by construction, Y,„ is the uniform 
amplitude of the raw M/L profile for x < Xi„, we 
expect it to be related to the "central" mass-to-light 
ratio of the galaxy, where by "central" is implied 
the distance at which the measurement of the cen- 
tral velocity dispersion (ctq) is obtained (at x = xq, 
with Xq typically less than x,„). However, in a real 
system, the true central M/L cannot be securely de- 
termined from the measurement of ctq alone, ow- 
ing to uncertainties about the validity of the assump- 
tions that are invoked, in order to translate knowl- 
edge of (To to that of mass enclosed within xq (us- 
ing virial theorem). Such uncertainties basically stem 
from the presence of anisotropy in phase-space. In 
the mass modelling trick advanced in Paper I, room 
is allowed for the accommodation of such uncertain- 
ties, as long as the deviations from the assumptions 
used in the virial estimate of mass are not atypically 
more than what has been observed with real ellipti- 
cals (Padmanabhan et a/.., 2004,) ; this is discussed in 
the following paragraph. 

In Paper I, the virial estimate of the central M/L, 
from CTQ was parametrised by a. It was found that 
for a measured ctq, as long as Y,„ lies within a range 
of values (the details of this range correspond to the 
given a), compatibility between the predicted and 
known (model) mass density distributions is ensured. 
In other words, for a calculated from a given ctq, Ti„ 
can be safely chosen to belong to a range of M/L 
values: Y™" to Y""". Such positioning of Y,„ can be 
checked by comparing the mass distributions recov- 
ered with Tin = T'"'" and Y,„ = T""", for consis- 

tn ni in ' 

tency. Here Y'"'" - a while Y""*^ is an unknown 

J m in 

function of the photometric parameters and a - say 
g(n,Xe,a). 



g(n, Xe, a). 



(3) 



where gin, Xg, a) is unknown and the choice of its 
dependence on the photometric parameters is the fol- 
lowing. In Paper I we had success upon choosing the 
functions / and g as dependent on the photometric 
details of the system; this was the case for a suite of 
model power-law galaxies. Such "success" is qual- 
ified in terms of the identification of the hypothe- 
sised dependence on the photometric properties of 
the model galaxies. Motivated by this, we endeavour 
to find forms of / and g. 

Thus, there are two unknown functions that we 
wish to constrain: /(«, X^, Y,„) and g(n, Xc,a). These 
functions, when known, will provide Tout from Y,„ 
which will be known from a, (i.e. ctq) and the ob- 
served brightness distribution of the galaxy. Once 
we know Tom from Y,„, we would then have fully 



4 



Dalia Chakrabarty and Brendan Jackson: Mass Distributions Only from Photometry & ctq 



characterised the raw MfL distribution over x. The 
smoothed out version of this raw MjL profile will 
then be advanced as the true MjL profile of the sys- 
tem at hand. This in conjunction with the luminosity 
density distribution will allow knowledge of the mass 
density distribution. 

We hope to recover the analytical forms of these 
two unknown functions for Sersic galaxies to 3Xe, 
through an analysis of a sample of model Sersic 
galaxy surface brightness profiles. Analytical fits are 
sought to the data that is collated from the suite of 
models that we work with, in order to recognise pat- 
terns, if any, that may show up in the relations be- 
tween the various quantities, in particular, Tour'^in 
and T,„ - a. The variation in these relations with 
changing models is then explored to unravel the re- 
liance of these relations on the photometric parame- 
ters. 



3. Models 

We identify the relation between Y,„ and To„, in 
Sersic galaxies, by recognising (and then quantify- 
ing) the patterns that emerge between these quanti- 
ties, as diverse model galaxies are scanned. To rem- 
inisce, Sersic galaxies are those, the surface bright- 
ness profile I{R) of which can be approximated as: 



/(/?) = /oexp[-/7„(W)'/«], 



(4) 



where /q is the central surface brightness, is the 
projected effective half-light radius, n is the Sersic 
index that determines the curvature of the bright- 
ness profile and bn is a functi on of n: for n > 
0.5, b„ ^2n-l/3 + 0.00987 6M dPrugniel & SimienI 
ll997l:lLima Neto et al.\l99% . Thus, the Sersic model 
represents a 3-parameter family. The correspond- 
ing deprojected spherica l luminosity dens i ty dis - 
tribution is discuss ed in [Terzic & Graham I ()2005l); 



iMazure & Capelatol (l2002|). In 
20071), he spherical models of 



Terzic & Sprague 



Terzic & GrahamI 



(|2005) were generalised to include triaxiality. The 
luminosity density profile in the triaxial case is ap- 
proximately given as: 



-b„(s/R,Y' 



where 



(5) 



0.6097 0.05563 
p = 1.0 + 5 — forne (0.6,10). 



Here, the equation for p is due to iLima Neto et al] 
(1 19991) and the ellipsoidal coordinate s is defined as 



n2 - 



(6) 



for the axial ratios of a:b:c for the ellipsoidal system 
at hand. We use in place of (see Section 7.5). It 
is to be noted that Equation|5]is an approximation an d 
other forms have been used bv lTruiillo et al.l (l2002h . 



As in Paper I, in our models, we assume 
oblateness and an inclination of 90°, with a uni- 
form projected axial ratio of 0.7; the work of 
iPadilla & StraussI (120081) corroborates such consider- 
ations of geometry and ellipticity. Also, the model 
systems are assumed to be viewed in the ACS z- 
band, at a distance of 17 Mpc. Though we recover 
the sought functions for these chosen configurations, 
generalisations to other systems will be suggested in 
Section |53]and Section|5.4| 

Actuallv lTerzic & Spraguj (l2007l) gives the more 
general 5-parameter mass density model which de- 
scribes Sersic galaxies with a core inside a given 
break radius. However, we work with the simpler 
case of Sersic galaxies that can be qualified by a 3- 
parameter model. In fact, we constrict this further, 
by normalising the central luminosity density po to 
1000 MqPc^^. Also, we look for the luminosity dis- 
tribution along the major axis. Thus, our models are 
distinguished only by the Sersic index « and the Xg 
defined in equation [T] 

Models were formulated for Sersic indices n=3, 
4, 5.3, 6, 7, 8 and 9 with X, of 726 pc, 1000 pc, 1400 
pc, 2000 pc. While these models were extensively 
explored to achieve the sought functions, we inves- 
tigated selected models with n lying in this range 
and higher Xg values. In particular, we are inter- 
ested in models with Sersic index <4, with larger 
Xg. Our viable models include galaxies with «=3, 
Xe <4 kpc and n-4, X^ <5 kpc. Such an extended 
size range holds special interest for the greater frac- 
tion of observed ellipticals that correspond to this 
low n. Additionally, models with other n were also 
found to work toX^- 10 kpc. This choice of models 
is supported by the results that 

- in the nearby Vi rgo cluster, the ACSVCS survey 
dFerrarese et al .l 2006) reports that nearly 90% 
of the targetted early-type systems (Sersic and 
cored-Sersic included) fall in the range of R^ <5 
kpc in the z-band (which is the waveband di- 
rectly comparable to our models). All the pro- 
gramme galaxies with n €[3,4] were found to 
have Re < 2.2 kpc while the single galaxy with 
n 6(4,5.3] and R^ >5 kpc, is not a Sersic galaxy 
but is a cored Sersic system. Here, we remind 
ourselves that the reported effective radius R^ is 
really the geometric mean of the extent along the 
photometric semi-axes. Thus, the extent along 
the semi-major axis is greater (by about a fac- 
tor of about 1.2, for an axial ratio of 0.7) than the 
reported Re- In other words, the surrogate for X^ 
is about 1 .2 times the values of R^ quoted earlier 
in this paragraph. Even when this factor in taken 
into account, the range of our models covers the 
programme galaxies of the ACSVCS. 

- measurem ent of effective radius is waveband de- 
pendent (.Temi et all l2008t iKo & Iml l2005h . so 
that for n=3 and 4, our limits on the model Xe val- 
ues in the z-band, are compatible with observa- 
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Fig. 2. Righit panel: trends in To,,, with T,„, with changing Sersic index n, from model galaxies with Xf=726 pc. The 
colour coding used for the different n values is: yellow for n=3, magenta for n=4, black for n=5.3, red for n=6, blue for 
n=7, green for ;!=8 and cyan for n=9. In the left panel, T™'" is plotted against a in broken lines; this plot is universally 
given by T""" = a. T""' is plotted against a for the different model galaxies used for the figure on the right, with the 
same colour coding. For a model at hand, at a given a, T;„ can be chosen from the range defined by the values of T™" 
and T'"'" . 



tions dFerrarese et alJl2006t rRuiillo et alJlIOOll: 
iLa Barbera et al.ll2005l)" 
- at high redshifts, systems display evolution to - 
wards higher compactness (Buitr ago et al.l2008l) . 
Though this evolution i s marked for s ysteni s 
at z >1.7 ( s ample of iBuitrago et all I2OO8I) . 
iFerreras et alj ( 1200 9*) suggest a milder size evo- 
lution for the most massive galaxies for 0.4< 
z <1.2. Thus, in the higher redshift range that 
we propose our scheme to be most useful in, if 
anything, the size range of our models would be 
more suitable. 

This corroborates our choice of models, particularly 
in regard to the values of Xg. 

Sersic indices 2 or lower appeared not to be vi- 
able for the formalism to function, implying that this 
formalism is suitable only for elliptical systems (see 
Section[3ll. Sersic indices greater than 9 were not ex- 
amined since such systems are very rare. 

These luminosity density models, described by 
Equation |5] and the used values of n and Xg, were 
embedded in an NFW-type dark halo jNavarro et al.\ 
Il996h of mass density pdark, to give a total mass den- 
sity of: 



Pr(x) = Pdark(x) + api(x), 

where by choice. 



Pdark(x) - 



4jtx{x + Ts) 



2 ■ 



(7) 



(8) 



Here Ms and are mass and length scales of the 
halo, respectively. Our results are valid for halo pa- 
rameters that correspond to the points in the green 
quadrilateral in the Ms - rs space that is depicted in 
the right panel of Figure [TOl 

We ascribe +l-cr errors of about 10% to the lu- 
minosity density distributions that we generate and 
search for compatibility between the model mass 
density distribution and the predicted one, within 
these error bars. 



4. Method 

Combinations of the parameters a, Tj,, and To,,, that 
showed compatability with the known (model) mass 
density profiles, up to 3Xe, at the aforementioned 
four separate (Ms, r^) coordinates, were searched for 
by our smoothing formalism that is automated. 

We record the list of T„, values that imply com- 
patibility for a chosen a (chosen typically in the 
range of 2 to 10). The starting value of T,„ was 
typically about a while the upper value was set as 
3a, which was always sufficient to find the whole 
range of T„, values that correspond to compatability. 
Similarly, we record the To,,, value corresponding to 
agivenT,-,,. 

The plots of Y""" and T""", as functions of a are 

w m ' 

monitored, with the aim of recognising the analytical 
form of these functional dependences. Similarly, the 
plot of To,,, as a function of T,„ is analysed at differ- 
ent n and Xe, to identify /(«, X,,, Ti„). 



6 



Dalia Chakrabarty and Brendan Jackson: Mass Distributions Only from Photometry & ctq 



5. Results 

With Xe held at 726pc, as n is varied, the plots of T,„ 
against a are shown in the left panel of Figure 2 while 
the relation between Tout against Y,„ is depicted for 
these models, in the right panel of this figure. 

5.1. Relation between T,„ & a 

When the data collated from the diff'erent models 
is plotted it is noted that while T""" = a, the 

^ in ' 

log(T™'") - log(Q') relation is well fit by a quadratic 
function of log(Q;). This latter functional form has 
only a weak bearing on the Sersic index n, though 
T™-* is found to increase slightly as n is increased 
from 5.3 to 9. Thus, by denying this increased ampli- 
tude of log T™'" with higher values of «, we merely 
constrict the range of values from which T,„ can 
be chosen. Thus, our attempt at simplification of 
the sought functional form preempts a small reduc- 
tion in the applicability of our formalism to more 
anisotropic galaxies than what is effectively allowed. 

In Figure |3] we present the values of log(T™') 
at distinct values of log(Q'), for different Xe, and a 
single Sersic index of 5.3. The functional form of 
the relation is recovered from our analysis, and this 
is over-plotted on the data, in solid lines. This func- 
tional form is given by: 

T™^ = Ao(X,) + A 1 (X,) log a + A2(X,)(log a)\9) 
Ao(Xe) = -0.00747849 H- 0.000360389X, 

-1.31416 X 10"X? 
AiiX,) = 0.597653 - 2.61661 x lO^^X, 

+9.34449 X IQ-^X^ 
A2{X,) = 0.396038 - 7.77987 X lO^^X, 

-3.03313 X lO^'^X;. 

5.2. Relation between T^,,, & T,„ 

The data from the assorted models indicate that Tout 
is noted to fall exponentially with Y,„ ; in fact, a good 
fit to this relation is given by the following equation: 

T„„, = A(n,X,)+Ao(n,X,)cxp[-rjT(n,X,)], (10) 

where A, Aq and r are functions of X,, and n. The re- 
lations between n and A, Aq and t are shown in the 
left, middle and right panels of Figure H) for the four 
different X,, values that we use. From our models, we 
seek the functional form of these relations; these re- 
covered functions are over-plotted in solid lines, on 
the data in the three panels in Figure|4] This compar- 
ison indicates that the known trends in log A, InAo 
and r(n), with increasing log(n) and n, for the as- 
sorted Xe values, are well replicated by the predicted 
functional forms of the relevant quantities. 




log(a) 



Fig. 3. The values of log T'."°' obtained at distinct 
values of log a, from our models, for «=5.3 and dif- 
ferent X/, in black filled circles for Xg-726 pc, red 
for Xe=1000 pc, green for ^^=1400 pc and blue for 
Xp=2000 pc. The predicted dependence of logT™" 
on log a (Eqn. |9]l is plotted in solid lines, for the 4 
different values of Xe that we use, in corresponding 
colours. This functional form holds approximately 
for other values of n as well (n > 3); the amplitude of 
YHMA Qjjjy increases slightly with n. 

We advance the following formulae for A, Aq and 

t: 

log[A(n,X,)] = Co(X,) + Ci(X,)log(n) where (11) 
CoiXe) = 0.18832 - 8.40018 X lO^'^X, - 

1.19449 X lO^X^ and 
Ci(Xe) = 1.15563 -H0.00207Z,- 
2.06978 X IQ-^X^ 

MAo{n,Xe)] = DoiXe) + Di{Xe)n + DiiXeW (12) 

+Di(Xe)n^ where 
DoiXe) = 480.078- 1.32685X, 

H-0.00106312X,2 

-2.63527 X IQ-'^xl 
Di(Xe) = -192.617 -t-0.537814X, 

-0.00043 1932X,2 

+ 1.07350 X 10"X^ 
D2(Xe) = 24.9361 -0.0699530X, 

+5.63996 X lO^X^ 

-1.40766 X lO^'^X,^ and 
D3(Xe) = -1.03518 + 0.00293420X,- 

2.37672 X IQ-^X^ 
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log(n) Sersic Index (n) Sersic Index (n) 

Fig. 4. Trends in the three different parameters that define the exponential fall-off^ of To,,, with T„,, with 
changes in n and X^; marked in filled circles of colours that are distinctive of the Xg value used - black for 
Xe-126 pc, red for Xe-IQQQ pc, green for Xc=1400 pc and blue for Ze=2000 pc. These trends have been 
spotted during the parametric analysis that we undertake with our model galaxies. The recovered functional 
dependence that defines these trends are marked in solid lines of colour corresponding to the value of X^,. 
These functional dependences were recovered using model galaxies with {« - 3, 5.3, 6, 7, 8, 9allXe < 2000pc) 



+5.96666 X 10 



-''xl 



T(n,X,) 



ToiXe) 



tq(X,) + Ti(X,)n + T2(X,)n^ + (13; 
TT,(Xe)rr' + T4(Xp)n'* where 
-3650.81 + 8.99277X, - 0.00649587^,2 
+ 1.43425 X lO^'^X^l 



1903.74 - 4.64150X, + 0.00331538X^ 
-7.22203 X 10"X^ 
T3(Xe) = -367.254 + 0.885327X,-0.000623864X, 

and 



needs to be scaled. In fact, the value of T,,,,, corre- 
sponding to observations in w should be scaled by 
I factor of lo^^^^^"^^-'. Here Lq^. is the solar ab- 
olute magnitude in the w-band. In the ACS z-band 
L0;=4.52. However, this transformation into the w- 
band is really an approximation since we assume all 
the way that Xg is th e same over wavebands. This is 
not true jTemi et al4.2008;.Ko & Im.2005.) . 



+ 1.33667 X 10"^X„3 



T4iXe) = 31.1895 -0.0742165X, + 5.14858 X 10" 



-1.08171 X IQ-^xl 



These relations are valid for all model systems 
with n >3 and that too, for n=3, only systems with 
Xe <4000pc and for n=4, X, <5000pc. 

5.3. Photometry in other wave-bands 

It is always possible that the observed photometry 
is presented in a wave-band different from the ACS 
z-band for which we predict the above relations be- 
tween the properties of the raw M/L, though the z- 
band is in general a better choice for high redshift 
systems than a bluer waveband. If the available pho- 
tometry is in a band different from z, (say the w- 
band), then the factor by which every value of T,„ 
and Tout should be changed, is obtained from the fol- 
lowing considerations. We realise that if w is such 
that a model galaxy is brighter in z than in the w- 
band, then the M/ L ratio by which the inner and outer 
parts of the brightness profile in z need to be scaled, 
are smaller than the same by which the profile in w 



5.4. Distance to the gaiaxy 

It is important to enquire about the stability of the 
posited forms of To,,, and T,„ when the distance to a 
Xgalaxy is different from what has been used in the 
models, namely 17 Mpc (approximate distance to 
Virgo). The only influence of the distance D (in Mpc) 
to the observed galaxy is in affecting the luminos- 
ity density distribution that is obtained be deproject- 
ing the observed surface brightness profile, through a 
term that is linear in D"' . It is this luminosity density 
profile that is scaled by T,„ and To„, in the inner and 
outer parts of the galaxy. Thus, if the galaxy is in re- 
ality further than 17 Mpc, the implementation of the 
suggested M/L values would amount to an overesti- 
mation of the luminosity density. To compensate for 
this, we need to modulate the relevant To,,; and T,„ 
values by the factor D/ll. 



6. Tests 

In this section we discuss the testing of our advanced 
parametric forms of To,,, - f{n,Xe,Ti„) (Equation 
9 with inputs from 10, 11 and 12) and T^^ = 
g(n,Xe,a) (Equation 8). The scheme is tested on 
models that were excluded from the fitting exercise 
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that resulted in our identification of the forms of the 
functions / and g. 

Figure |5] shows a comparison of the predicted 
and model mass density profiles along the x-axis for 
the model galaxy that has a Sersic index of 3.5 and 
X(,=3000 pc. a is set to 5 => ctq a;250 kms"', av- 
eraged over a radius of about X(,/8=375 pc. For this 
system, we extract T™" a^S.S and Yo„, ^0.5. This 
pair of Tout - Tj,, values are used to scale the luminos- 
ity density distribution for this model galaxy which 
is then smoothed to offer a final mass density distri- 
bution along the major axis. These To,,, - T„, values 
imply a predicted mass density distribution (in red in 
Figure |5j that tallies favourably with the model (in 
green in Figure |5]l. 

The predicted mass density distribution when 
compared to the luminosity density distribution of- 
fers the final or smoothed M/L which is shown in 
Figure |6] to be significantly different from the raw 
M/L profile (in this black line). 

However, this was the example of one given 
model galaxy, characterised by a given central ve- 
locity dispersion. Tests were undertaken to validate 
the relations predicted between the properties of the 
raw M/L and the photometric parameters of model 
galaxies, across the full range of ctq. The tests were 
carried out with model galaxies with Sersic index 
of 4 and X, of 726 pc and 1400 pc. In Figure [T] 
we see the calculated Tout value corresponding to a 
given T,„, for which compatibility is noted between 
the known (model) and predicted total mass density 
distributions, to 3Xe. Such To„, values are shown in 
the two panels of this figure, in black dots. 

It is to be noted that these calculated values of 
Tout very closely straddle (within errors of ±5%) the 
analytical relationship between Yo„, and T,„ that is 
predicted in Equations 9, 10, 11 and 12, for given Xe 
and n; in this case for «=4, Xe-126 pc (left) and n-A, 
Xe=1400pc (right). Thus, these tests offer confidence 
in the formalism that we suggest. 



7. Applications 

In this section, we check out the efficacy of the ad- 
vanced scheme in recovering gravitational mass den- 
sity distribution of two elliptical galaxies, NGC 4494 
and NGC 3379, to 3Xg. Our predicted mass distribu- 
tions are compared to independent dynamical mass 
models for these systems. 

7.1. NGC 3379 




NGC 337 9 was reported by iRomanowskyefa/.l 
((2004) and 'Doug las et. all (l2007 h to contain very lit- 
tle dark matter on the basis of a Jeans equation 
analysis of the kinematics of around 200 plane- 
tary nebulae (PNe) that reside in the dark halo of 
this galaxy. An independent estimation of the dis- 
tribution of the total mass density of this galaxy 



10 100 
Distance along x (pc) 

Fig. 5. Figure to compare the predicted total mass 
density distribution with x (in red) with the known 
(model) mass density profile (in green), for the model 
described by n=3.5 and X^=3000 pc. For an a of 5, 
this model yields T.™^ ;*5.8 and Y^,,, ;*0.5. The pro- 
file in red results on scaling the luminosity profile of 
this model galaxy by the raw M/L distribution that 
is defined by these "inner" and "outer" amplitudes 
and subsequently smoothing the discontinuous mass 
density profile thus obtained, (using the smoothing 
prescription mentioned in the text). 



was performed by Chakrabarty, 2009 (submitted to 
Af), by implementing these PN e velocities in the 
Bayesian algorithm CHASSIS (Chakrabait v & Sahd 
[2001; Chakrabarty & Portegies Zwart 2005). As ac- 
knowledged by Chakrabarty (2009), these mass 
distributions from CHASSIS indic ate somewhat 
higher masses than the estimates of iDouglas et. all 
(12007 1. owing to the assumption of isotropy within 
CHASSIS. While details of such mass estimation 
techniques are irrelevant to the current work, here 
we present a comparison between the mass density 
profile obtained from our formalism with the same 
obtained from CHASSIS. We also present a compar- 
ison between the cumulative mass result M{r) via the 
quantity defined as Vc - ■\/GM{r)/r, where vv is re- 
ferred to as the circular velocity. 

Our predictions are based on the analysis of 
the photometry of NGC 3379 in the B-band 
dCapaccioli et all 1 19901) and the data for pro- 
jected central velocity disp ersion, as given by 
'Statl er & Smecker-Han3 (119990 . Using this o-p data 
along the major axis of the galaxy, (Table lA of 
IStatler & Smecker-Hand 19991) . we get an a of about 
7.5. The photometry suggests an X^ of 2.2 kpc. 
Additionally, we deproject the B-band surface bright- 
ness profile using the Bayesian deprojection algo- 
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Fig. 7. Comparison of model To„,, as obtained from our calculations (black filed circles) and predicted To,,, 
(red line), for given Y,„, corresponding to model systems with Sersic index = 4 and Xg-126 pc (left) and 
Xe-\AQQ pc (right). The predicted curves are obtained by using the values of n and X^, for the model test 
system, in Equations 9, 10, 1 1 & 12. The calculated T„„, values result from our analysis. 
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Fig. 8. The radial distribution of the recovered total mass density is depicted in black, on the left. On this 
is superimposed the mass density profile obtained by implementing the kinematic information of a sample 
of planetary nebulae in NGC 3379 (the same sample as used by Douglas et. al, 2007), in the Bayesian 
algorithm CHASSIS (Chakrabarty 2009, communicated to AJ). The quantity Vc, as recover ed from these 
mass distributions are shown in corresponding colours on the right, with the Vc profile from lDouglas et. al 
(120071) overplotted in red. 



rithm DOPING jChakrabartv & Ferraesell2008l) . We 
consi der NGC 3379 to be at a distance of II kpc, 
as in iDouglas et. all ( l2007h . Our predicted values of 
(upper limit on) T,„ and Tout, as modulated by differ- 
ences in the wave-band of the available photometry 
and distance to the system, are about 9.3 and 5.9. 

The resulting mass density distribution that we 
advance for NGC 3379 within 3X(., is depicted in 
black on the left of Figure[8] This is compared to the 



mass model identified by CHASSIS. Ve derived from 
our calculated mass distribution, from the dynami- 
cal mass modelling by CHASSIS and that reported 
by Douglas et. al ( 2007.) are represented in the right 
panel of Figure [8] 

The comparison of the total mass density pro- 
file indicates the clear trend for our predicted 
mass distribution to be on the lower side com- 
pared to the gravitational matter density provided 
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Fig. 6. Figure to compare the final total (local) MjL 
distribution with the initially chosen raw Mj L profile 
(in broken lines), for the test galaxy the total mass of 
which is presented in Figure|5] The final MjL is pre- 
sented within 3Xe. The raw MjL profile is the same 
as that shown in Figure [T] as is apparent from this 
figure, it is significantly different from the final MjL 
profile. 



by CHASSIS. This is only to be expected since 
CHASSIS in its current form, assumes isotropy, 
which affects CHASSIS idiosyncratically to spuri- 
ously enhance mass density, as acknowledged by 
Chakrabarty (2009). Consequently, the profile ad- 
vanced by our work is also on the lower side of 
the profile that follows from CHASSIS, though our 
result compares bett er with the result advanced by 
iDouglas et. ail(l2007l) . 



7.2. NGC 4494 

NGC 4494 is a nearby elliptical, the mass dis- 
tribution of which to 7 effective radii has re - 
cently been presented by iNapolitano et alj ( l2008l) . 
Th e distance to this galaxy is given as 15.8 Mpc 
by INapolitano et all feoOSi) . The V-band photom- 
etry of this galaxy is presented in Table Al of 
INapolitano e t al. (2008). This surface brightness pro- 
file indicates an X^, of 3.37 kpc and is deprojected, 
given the radial variation of the projected eccen- 
tricity (assuming an oblate geometry and edge-on 
viewing). The projection of such recovered luminos- 
ity density distribution is shown in the left panel of 
Figure |9] compare d to the surface br i ghtnes s data of 
this galaxy (in red) . INapolitano et aTl (l2008l) also cite 
the cen tral velocity disper sion of this system as 150 
kms~' dPaturel et al.l2003T) . Considering this cr,, to be 



the dispersion averaged over Re/S, where Re-48".2 
for this galaxy, we obtain an a of 2. Our predicted 
values of T,„ and To„, are about 2 and 0.121. These 
are used to characterise the raw M/L, which when 
smoothed, gives rise to the total mass distribution 
that is shown in the middle panel of Figure |9l 

The total M/L ratio of NGC 4494, in the V-band 
is shown on the right panel of Figure |9] while the 
enclosed gravitational mass distribution is depicted 
in the middle panel. This is very similar to the ra- 
dial mass distribution of NGC 4494, within t he in- 
ner 3Xe, as reported bv lNapolitano et al] (l2008l) (their 
Figure 13). 



8. Discussions and Summary 

We have presented a novel mechanism for estimating 
total mass density profiles of elliptical galaxies that 
can be described by a Sersic-type surface brightness 
distribution. This formalism uses nothing in excess 
of what is typically available in the observational do- 
main, namely, photometry and the central velocity 
dispersion profile. This allows the implementation of 
this scheme even for systems at high redshifts. (Such 
an implementation of this scheme will be presented 
in a future contribution - Chakrabarty & Conselice, 
under preparation). The advanced scheme uses the 
methodology presented in Paper I, in the context of 
Sersic galaxies. 

To begin with, we generate a sample of luminos- 
ity density distributions that project to Sersic bright- 
ness profiles, of assorted values of n and X,,. These 
toy galaxies are then assigned various values of the 
local central mass-to-light ratio (parametrised by a). 
For each such configuration, we monitor the allowed 
range of T,„ and the value of T,,,,, that corresponds to 
any such Y,„. This pair of T,„ - To„, values define the 
raw or unsmoothed two-stepped local mass-to-light 
ratio distribution with x. The local mass density dis- 
tribution in the system (to 3Xg) is judged by scaling 
the luminosity density profile by the raw M/L pro- 
file and then smoothing the result by the prescribed 
smoothing routine (two successive applications of a 
box filter of size corresponding to X^). 

The dependence of To„, on T,„ is found to be well 
approximated with an exponential decline that is de- 
scribed by three free parameters, which are actually 
all functions of the model characteristics, namely, n 
and Xe. The dependence of these three parameters on 
n and X^ are extracted from our model-based study 
and analytical fits to these trends are presented. Such 
fits are noted to be favourably represented by poly- 
nomials. 

The reliance of the band of the allowed values of 
Y,„ on a is also quantified as ranging from a to T™^, 
where the form of the function TJ.™^(n, X^, a) is also 
found to be polynomial in nature. 
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Fig. 9. Figure to bring out the results for the galaxy NGC 4494. The left panel depicts a projection of the 
luminosity density distribution that we estimate for the galaxy (in black), compared to the V-band surface 
brightness data (in red). The total gravitational mass distribution of the galaxy is represented in the middle 
panel while the mass-to-light ratio on the V-band is shown on the right. 



8.1. Deviation from predicted trends? 

Even though the polynomial fits are found to do a 
good job for most configurations, for certain mod- 
els, the predictions appear incorrect (by no more than 
10%), eg. T may be judged under-estimated for «=8 
and 9 for the configuration: X(,=2000 pc. However, 
this should be interpreted as erroneous recording of 
Tout for these T,„ values from our models, rather 
than a failure of our predictions. Typically, such er- 
rors emanate from lack of refinement in the step-sizes 
used in our work and is potentially amendable. 

8.2. Applicable to. . . 

As delineated in Section[3] the formalism that we ad- 
vance here, is applicable only to galaxies with Sersic 
indices n > 3. It is only for such systems that the 
luminosity density profile is steep enough, i.e. falls 
quickly enough to ensure that in the "outer" parts 
of the system the luminous matter density is com- 
pletely overwhelmed by the dark matter density, so 
that the total matter density is effectively a represen- 
tation of the dark matter distribution. Thus, changes 
in the central M/Lor a do not affect the M/L in these 
"outer" parts. Now, to qualify "outer", we state that 
this region corresponds to x ~ 3Xc. This is anyway 
true, only if the luminous matter fraction is not too 
small, in which case, we need to settle for progres- 
sively lower M/L in the outer parts, as a increases. 

8.2.1. Effect of lowering Sersic Index 

If however, the luminous density profile is too flat, 
(as for n < 3), the contribution of the luminous mat- 
ter, to the total mass, is significant, even in "outer" 
parts. In order to achieve similarly shaped total mass 
profiles, the value of To,,, will need to increase with 
increasing a. Then, the exponential fall-off of To,,, 



with T,„ will not be noted. Qualitatively speaking, 
this is the reason why our formalism will not be valid 
for galaxies with Sersic indices less than 3. 

8.2.2. Effect of increasing size 

However, this increased flattening of the luminosity 
density profile can be attained even for n >3, for 
large values. Thus, in such systems, for the same 
reason as in the last paragraph. To,,, will increase 
with an increase in the inner M/L instead of falling 
exponentially with increasing T,„. Thus, based on 
these considerations, on the basis of our experiments, 
we find that the following galaxies cannot be accom- 
modated within our suggested formalism: 

- n=3 and X, >4000 pc, 

- n=4 and X, >5000 pc, 

- n <3. 

8.2.3. What the ranges mean 

However, it is important to keep in mind that these 
ranges have been ascertained on the basis of a dis- 
cretised scanning of the parameters n and Xe, of our 
two-parameter galaxy models. Thus, when we say 
that models with n <3 cannot be explained with 
the predicted relationships between a, , T,„, To,,?, 
what we really mean is that our experiments with 
a model with n-2 indicated failure of adherence to 
these predicted trends. It is very much possible that 
a model with an intermediate Sersic index can be 
explained by our predictions, such as a galaxy with 
n=2.3, say. In other words, the ranges of n and Xg de- 
scribed above owes, to our dealing only with models 
ascribed the aforementioned discrete values (values 
mentioned in Section |3]l. Consequently, the ranges 
listed above are a conservative interpretation. 
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The practical question to ask is, are these sug- 
gested ranges compatible with real galaxies? As 
itemised in Section [3] our choice is supported by 
the trends observed for the Virgo cluster, within 
the ACS Virgo Cluster Survey, in alliance with 
the predicted evolution towar ds greater compact - 
ness in high-redshift systems ( Ferreras et alJl2069h . 
Furthermore, keeping in mind that a waveband de- 
pendence of effective radius exists, our model val- 
ues in the z-band faU within lirnits indicated by ob- 
servations (Ferrarese et al.ll2006l : lTruiillo et al.ll200ll : 



iLa Barbera et alj|2005l) 



Additionally, we question the relevance of the 
upper limits on the effective radii at «=3 and 4, 
which are smaller than the range permitted for higher 
values of n. As far as the correlation between the 
shape parameter n and log of effective radius is con- 
cerned, it is well known that in g eneral, half-light ra- 
dius decreases w ith decreasing n (iBoselli et al.ll2008c 
Cote et al.ll2008HNaab & Truiilloll2006l:lBrown et all 



2003l:lTruiilloetalJl2001h . In fact, such a correlation 



is physical and is hinted at by the global relations 
between structural parameters; it is not merely a re - 
flection of the fitting procedure (ITrujillo et alJl200ll) . 



8.2.4. Effect of model dark matter distribution 

Of course, it could be perceived that such fate of the 
mass distribution at x ~ 3Xp would be dictated by 
the exact details of the underlying dark matter distri- 
bution that defines the model galaxy. Indeed, for the 
high values of «, the M/Lat x > 3Xf (T,,,,,) is dictated 
by the details of the dark matter distribution, but it 
is the smoothing of this raw M/L profile that gives 
us the end product, namely the total local M/L ratio 
profile of the galaxy upto the benchmark x of 3Xc. It 
is precisely this smoothing procedure that brings into 
the final M/L structure, elements of the mass distri- 
bution in the inner 3Xc in the galaxy, as well as a 
relatively less significant signature of the same that 
Ues outside the radius upto which the distribution is 
sought. The difference that smoothing makes to the 
raw M/L profile, as distinguished from the final M/L 
profile, is clear in Figure|6] 

Thus, given such (1) implementation of our 
smoothing procedure and that (2) the jump radius in 
the raw M/L profile is itself 3Xe, we find that our ex- 
periments bear the fact that changing the mass scale 
and length scale of the used NFW dark matter density 
distribution across the wide ranges do not affect the 
mass configuration for x < 3Xc, though further out in 
the galaxy, the influence of changing the DM density 
distribution picks up quickly. The consistency in the 
Tout - '^in relation, noted with changes in Mj and rj 
is brought out in the left panel of Figure [TO] The al- 
lowed ranges in and are provided in the right 
panel of Figure [TOl 



8.2.5. Wliy clioose NFW? 

The modelling of the dark matter distribution in early 
type galaxies is obviously indirect and consequently 
difficult and unreliable. Ferreras et al. (200 8) suggest 
a large scatter in the outer slope of the dark mat- 
ter distri butions that they reco ver for their targeted 
galaxies. iGavazzi et alj ( l2007l) corroborate an NFW 
modelling of the dark haloes of their sample galax- 
ies, suggesting that the total matter density tends to 
an overall isothermal form. Given this degree of un- 
certainty, it might be argued that there is not much 
sense in splitting hair to decide between an NFW and 
isothermal dark matter distributions, as long as the 
dark haloes that we probe in our analysis are com- 
patible with observations or simulations. 

This is indeed the case, when we compare 
our ranges of halo characterist ics (right panel of 
Figur e [TOl ) to the suggestion by iKleinheinrich et all 
(12006 ) that the recovered range of virial masses of 
their NFW model to be [3.9 x 10", 7.1 x 10"] Mq. 
Such a mass range is compatible with the sugges- 
tion bv lHoekstraetal] ( 12004 ) for a fiduciary galaxy 
of luminosity Lb - 10^° L^q that has an NFW pro- 
file (M200 e ri0,6.4] X 10"M h, as indicated by 
Figure 4 in lHoekstraetani2004 . This mass range 
corresponds to a range of about 8 kpc to about 30 
kpc for r,. In fact, we cover these ranges in the runs, 
the results of which are presented in Figure [TOl and 
also scan haloes of lower masses. For various haloes 
defined within the red and black quadrilaterals in 
Figure [TOl compatible T,,,,, - T,„ relations were re- 
covered. 

This explains our choice of sticking with the 
NFW prescription for our dark halo model. For rea- 
sons similar to what we explain above, we suggest 
that upon smoothing, haloes of varying shapes but 
similar masses within 3Xe are not expected to affect 
the mass distribution within 3Xf, in the massive ellip- 
ticals that we deal with in this formalism. 

Thus, we see that our methodology is valid, irre- 
spective of selecting haloes defined by widely dif- 
ferent Ms and r, values, (i.e. NFW haloes defined 
by points inside the black and red quadrilaterals in 
Figure [Toll. This motivates us to choose to work with 
a smaller subset of all the haloes explored in these 
runs; in particular, we choose to work with haloes de- 
fined by the more massive half of the log(Mj) range 
that limits the red quadrilateral. This chosen subset in 
the Ms - rs space is bound by the green quadrilateral 
in Figure [TOl 

8.2.6. Effect of model geometry 

In addition to details of the DM distribution, we ad- 
mit that the model structure does include the free pa- 
rameters that describe the Sersic model of density 
distributions in general triaxial galaxies, as it ide- 
ally should, though this work has full potential of 
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Fig. 10. The left panel depicts the relation between To„, and T,„, obtained using the model with n-5.3 
and Ze=1000pc, from two different runs that seek this relation, given NFW haloes characterised by (M,, r,) 
values within the quadrilaterals in black (see right panel) and in red (see right panel). When the halo density 
parameters are as contained within the black quadrilateral, the obtained To,,, - T,„ relation is shown in open 
black circles in the left panel. When the halo density is characterised by and rj values as points inside 
the red quadrilateral on the right, the Tout - Y,„ relation is shown in red filled circles in the left panel. The 
predicted To„, - T,„ relation in shown on the left in solid black lines. The green quadrilateral represents the 
subset in the Ms - space, which defines the dark matter density distributions for which concurrent / and g 
functions were recovered in general runs done with assorted models, within our scheme. 



being extended into two-dimensions and include an 
axisymmetric description of the galaxy. Now, the tri- 
axiality of the models needs to be defined in terms of 
a chosen geometry and intrinsic eccentricities - it is 
of course not possible to constrain such characteris- 
tics from observations alone. 

In our models, these were fixed as oblateness 
with an axial ratio of 0.7. It is indeed possible that 
the unknown functions / and g that we attempt to 
constrain, harbour dependence on these intrinsic ge- 
ometric factors. However, the preparatory assump- 
tions involved in the deprojection of an observed 
brightness distribution, namely the underpinning of 
intrinsic geometry and viewing angle, are essentially 
unconstrained, unless the system is favourably in- 
clined or flattened. In other words, the unknown in 
our modelling are the usual quantities that render de- 
projection non-unique. 

As for the specification of the ellipticity in 
our models, a ty pical value has been adopted - 
iPadilla & StraussI ( 1200 8) suggest an axial ratio dis- 
tribution for a large sample of SDSS ellipticals, with 
a mode in the range of 0.6 to 0.8. Inspired by this, we 
use an axial ratio of 0.7 in our models. Again, the un- 
certainty in intrinsic ellipticity, cannot be known for 
a general observed galaxy. The quantification of the 
deprojection effects on the formulae provided above 
is possible, at least in a statistical sense, and it is en- 
visaged that the same will be pursued in the future. 



8.3. What if central super massive black holes? 

The case of a central mass condensation in the sys- 
tem was dealt with in Paper I, in reference to the ex- 
ample of the galaxy M87 - the recovered mass dis- 
tribution of M87 was demonstrated to be consistent 
with the same obtained from kinematical consider- 
ations. If independent measurements indicate an ob- 
served system to harbour a massive central mass con- 
densation, then the scheme delineated in Paper I will 
be followed. 

8.4. Why X, instead of Re? 

The usage of Xe instead of Re is preferred since Xe as 
per its definition here, as well as in Paper I, is derived 
solely from the inputs to the methodology, namely 
the surface brightness profile along the semi-major 
axis which we consider to be along the x-axis. Thus, 
any changes to the shape of the brightness profile will 
be directly reflected in a linear change in Xe but not 
necessarily so in the half-light radius that is estimated 
from isophotal analysis. In fact, here we use the same 
definition for Xe as in Paper I, (see Section 2). 

The usage of Xe should not be cause for concern 
since Xe is merely one definition of the semi-major 
axis effective radius and reduces to the conventional 
definition of the major-axis effective radius for n=43- 
When an unknown galaxy is being analysed within 
this formalism, its half-light radius is just as much 
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an unknown as is the that we define here. Thus, 
there is no loss of connection with observations by 
the implementation of X^. 



8.5. Effect of choice of smoothing prescription 

It merits mention that other choices for the smooth- 
ing prescription and x,„ may also work, but here we 
concentrate on the above mentioned configuration 
and the specification of T,„ and Tout are accordingly 
unique to these choices. 

8.6. Effect of choice of estimation of central local 

MIL 

In a similar context, it may be argued that the defi- 
nition of a that we use herein will leave an imprint. 
We equate a to the central local MjL, inspired by 
this result that is achieved in Paperl. Here the total 
gravitational mass found enclosed within the radius 
xo is M(x()) where M(x()) is linked to the 3-D velocity 
dispersion at xq via: o-y - GM(xq)I xq. Then, accord- 
ing to our definition, a is given as the ratio of Mix^) 
and the enclosed light within xq. For other definitions 
of a, other ranges of Y,„ will be valid (and therefore 
other forms of dependences of To„, on T,„, for the 
same galaxy). Thus, the pairs of / and g functions 
that we advance here, work for the used choice of the 
definition. 

We would like to emphasise that the extraction of 
the exact value of the mass enclosed within xo ( and 
therefore of a), from a measurement of ctq, is not 
the point of this exercise; uncertainties in this extrac- 
tion do not undermine the advanced results either, as 
long as the galaxy at hand is "not too" aspherical or 
anisotropic at xo. Here we qualify "not too" as those 
configurations for which we obtain consistent mass 
profiles using the two extreme values of T„,, that are 
allowed for the extracted value of a. Thus, when the 
method fails, we know that it does. As long as this 
aforementioned consistency is noted, choosing T,„ 
from anywhere within the range corresponding to the 
given a will lead to consistent total mass distribu- 
tions within 3Xe. Additionally, this range is neither 
too constricted nor too relaxed, as was discussed in 
Paper I. 

8. 7. Comparatively better applicability 

The presented device is based upon conclusions that 
are drawn from a sample of model Sersic galaxies. 
This would naturally imply that the success of this 
formalism is crucially dependant on the generality 
of the models. In particular, we have discussed the 
ranges of n and X^ for which our advanced results are 

' Semi-majo r axis effective radius h as been used before, 
for example bv lNaab & Truiilld (l2006h 



true. We have also discussed the effect of changing 
properties of the dark halo that we use in the models 
and find the advanced scheme robust to such model 
parameter variations. 

In fact, the formalism presented above is unique 
in its scope and structure. Estimates of total mass 
distributions in distant elliptical systems are diffi- 
cult and therefore rare; the formulae presented herein 
are therefore advantageous and could be treated as 
guides to decipher the total mass distribution of 
Sersic galaxies in large surveys. 

Most importantly, the advanced methodology is 
successful within a severely constricted data domain, 
compared to any other scheme that aims to obtain 
mass distributions in elliptical galaxies. All that the 
advanced method demands in terms of data is what 
is typically available - surface brightness profile and 
a measure of central velocity dispersion. The unde- 
manding nature of our method renders it applicable 
even at high redshifts. The simplicity of implementa- 
tion of the input data is advantageous in that it allows 
for the scheme to be used in an automated way, to ob- 
tain mass distribution s for large samples of galaxies. 
Tricks such as this and lNipoti et al.l (l2008h exploit the 
basic configuration within galaxies and offer novel 
ways for characterisation of distant systems. 
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